Force dependent fragility in RNA hairpins 



O 
O 

(N 

C 

On 



M. Manosas 1 , D. Collin 2 and F. Ritort 1 

1 Departament de Fisica Fonamental, Facultat de Fisica, 
Universitat de Barcelona Diagonal 647, 08028 Barcelona, Spain 
2 Merck & Co. Inc., Automated Biotechnology Dpt., North Wales PA 19454, USA 

(Dated: February 6, 2008) 

We apply Kramers theory to investigate the dissociation of multiple bonds under mechanical 
force and interpret experimental results for the unfolding/refolding force distributions of an RNA 
hairpin pulled at different loading rates using laser tweezers. We identify two different kinetic 
regimes depending on the range of forces explored during the unfolding and refolding process. The 
present approach extends the range of validity of the two-states approximation by providing a 
theoretical framework to reconstruct free-energy landscapes and identify force-induced structural 
changes in molecular transition states using single molecule pulling experiments. The method should 
be applicable to RNA hairpins with multiple kinetic barriers. 

PACS numbers: 82.35.-x,82.37.-j,87.15.-v 
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Single molecule pulling experiments allow to exert me- 
chanical force on individual molecules such as nucleic 
acids, proteins and macromolecular complexes []J. By 
recording force-extension curves it is possible to deter- 
mine free energies and kinetic parameters of biomolecules 
and search for intermediates and pathways in biochemi- 
cal reactions. Over the past years single molecule tech- 
niques have been successfully applied to investigate the 
breakage of molecular bonds in many biological systems 
such as proteins |2|, DNA molecules |3j, RNA molecules 
0, ligand-receptor binding Q and beyond, e.g. metallic 
gold nanowires stretched with AFM [fjj. Under mechan- 
ical load all these structures yield at different values of 
the applied force in a dynamical process that is stochas- 
tic and loading rate dependent. The study of breakage 
forces under nonequilibrium conditions is known as dy- 
namic force spectroscopy 0. A detailed comprehension 
of the rupture kinetics of biomolecular complexes has im- 
plications in our understanding of their kinetic stability 
which is important in enzymatic and/or regulatory pro- 
cesses. 

Here we investigate the unfolding/refolding kinetics of 
RNA hairpins using laser tweezers 0, Q . The RNA se- 
quence and its native structure are shown in Fig^ To 
manipulate the RNA hairpin two beads are attached to 
the ends of the RNA hairpin by inserting two hybrid 
RNA/DNA handles One of the beads is immobilized 
on the tip of a micropipette while the other bead is cap- 
tured in the optical trap. By moving the micropipette a 
force is exerted upon the ends of the RNA hairpin and the 
force-extension curve (FEC) recorded. In Fig^we show 
the experimental FEC corresponding to a complete cycle 
of a ramping process where the force is first raised and 
relaxed afterwards. In the pulling process the molecule 
is initially in its native folded structure, and the force is 
increased at a certain rate r until the molecule unfolds. 
If the process is reversed, i.e. the force is decreased at 
rate —r, the molecule folds back again (relaxing process) 
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FIG. 1: Experimental FEC for the RNA hairpin in a ramping 
cycle where the force is first increased and after decreased at a 
rate of r = 3pN/s, the blue and green curves representing the 
pulling and relaxing processes respectively. The arrows indi- 
cate the unfolding and refolding forces, which correspond to 
the first unfolding and refolding events. The first part of the 
FEC corresponds to the stretching of the RNA/DNA hybrid 
handles. Inset: Experimental distribution of the unfolding 
(circles) refolding (triangles) forces. Number of pulls: 129, 
385 and 703 pulls for r=1.5, 7.5 and 20 pN/s (red, blue and 
green) respectively. 



The unfolding/refolding of the molecule can then be 
identified as force-extension jumps observed in the FEC. 
By repeatedly pulling the molecule many times we obtain 
the distribution of unfolding (refolding) forces, i.e. the 
force at which the first unfolding (refolding) event occurs 
along the pulling (relaxing) process (Fig^). The experi- 
mental distribution of the unfolding (u) and refolding (r) 
forces /9 u ( r ) (/) at different loading rates r is shown as an 
inset of FigQ] 

To model the hairpin we follow Cocco et al. ^lj and 
restrict the number of configurations of an N base-pair 
(bp) RNA hairpin to the set of configurations where the 
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first n bps are opened and the last N — n are closed (the 
total number of configurations being N + 1). The end- 
to-end distance is a well defined reaction coordinate for 
the unfolding/refolding reaction. We use the variable n 
to label the state of the hairpin; e.g the folded (F) state 
corresponds to n — and the unfolded (UF) state to 
n = N. The stability of each state n depends on its free 
energy, G n (f), at a given applied force / [111 ]. 



G n (f) = G°(n)+g(nJ) 



(1) 



where G°(n) is the free energy of formation at zero force 
and g(n, /) is the force-dependent contribution to the 
free energy. The latter is given by g(n, f) = G sa (n, f) — 
fx n (/) where the functions x n (/) and G ss (n, /) are the 
end-to-end distance and the entropic correction to the 
free energy of an M-bases long single stranded RNA (ss- 
RNA) at force /, M being the number of bases released 
after the opening of n bps. The latter can be computed 
as the reversible mechanical work needed to stretch the 
ends of an M-bases long ssRNA a distance x n (f), 



G ss (n,f) = 



*»(/) 



ssRNA 



{x)dx 



(2) 



where F ssRNA (x) is the FEC of the ssRNA .12\. The free 
energy landscape of a hairpin, G„(/), as a function of 
n is known to be rugged with different kinetic barriers 
(or transition states) depending on the sequence and on 
the applied force / We use the Mfold prediction 

(lif to extract the free energy of the molecule at zero 
force, G°(n). In what follows we take the F state as 
the reference state for the free energy, i.e. G n= o(f) — 
G°(n = 0) = 0. 

The kinetics of unfolding (i.e. the transition be- 
tween the completely folded (F) and the completely un- 
folded (UF) states) is an activated process with a force- 
dependent effective barrier, B c g (/), measured relative to 
the F state. The rates of unfolding and refolding, k u (f) 
and fc r (/), can be computed as the first passage rates 
for a Brownian particle to cross a force-dependent 
effective barrier B e g(f): 



k u (f) = k e-? B ««^-Mf) = k Q e~^ B ^~ G »^\ (3) 

where Gjyif) is the free energy difference between the 
F and UF states at force /, f3 — l/k B T with k B and T 
being respectively the Boltzmann constant and the bath 
temperature, and fco is an attempt frequency. An analyt- 
ical expression for B c g(f) can be derived from Kramers 
theory applied to the dissociation of consecutive bonds 
under mechanical force in the stationary approximation 
7. 16], 



4 T (/) = k B T\n 
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PGM) 



n=0 



(4) 



with h(n) — X)n'=o e_ ^ G ™ The variation in force of 
the effective barrier gives information about its position 
along the reaction coordinate: 



ggeffCf) 

df 



d(B cS (f)-G N (f)) 
df 



,(5) 



where x^ s (f) and x^g(f) are the distances from the ef- 
fective barrier to the F and UF states respectively. The 
location of the barrier along the reaction coordinate is 
related to the fragility of the molecule which determines 
how much the unfolding/refolding kinetics is sensitive to 
the force. To characterize the fragility we introduce the 
parameter fi defined as: 



x F cB (f)~^j(f) 



(0) 



[i < corresponds to a brittle structure -e.g. the case of 
hairpins stabilized by tertiary contacts where the barrier 
is located near to the F state-, whereas [i > represents 
a flexible or compliant structure, i.e. molecules that can 
easily deform under applied force and the barrier is close 
to the UF state [r?j . 

The rates are related to the force distributions by 



the expression p u (r)(/) 



)(/) 



exp[-<+) 



dy] 



with /' being the initial force in the pulling (relaxing) 
process 0|- The unfolding (refolding) rates read as 
K(r)(f) = Pu(r)(/) rPu( 1 r)(/) , where P u ( r )(/) is the proba- 
bility that the molecule remains in the F (UF) state along 
the pulling (relaxing) process until reaching the force /, 
Pu(r)(f) =+(-)/*/ Pu(r)(y)dy. Note that the experimental 
FECs show a force jump 5f (Fig^l when the molecule 
unfolds or refolds that corresponds to the relaxation of 
the bead in the trap after the sudden increase or decrease 
in the RNA extension. To compensate for this effect we 
shift the value of the folding forces by an amount equal 
to +(-)Sf/2. From (11 IMP and the unfolding (refolding) 
force distributions, p u ( T )(f), we can extract the effective 
barrier as: 



Pu(/> 



Pu(f)k 



(7) 



Pr(/)7 



(8) 



^ P (/ + ^)^°(iV) + ^,/)-Iln Lpr(/)fco 

Using polymer theory f3| we can estimate g(n, f) so 
the expressions I|7I8|I have only two unknown parame- 
ters, G°{N) and k . We determine G°(N) by collaps- 
ing into a single curve the effective barrier estimates 
and (JHl corresponding to the pulling and relaxing 
processes at different loading rates. From our data we 
get G°(N) — 64.5fcsT in very good agreement with the 
Mfold prediction, G^ fold (iV) = 38Kcal/mol = 63k B T 
[l5| . In Fig|3we show the force dependent effective bar- 
rier obtained using this method. We then determine 
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FIG. 2: Effective barrier B°g P (f) I0H} from the pulling 
(circles) and relaxing (triangles) experimental data compared 
with the prediction by Kramers theory, B^f (continuous 
line). Inset: Free energy landscape Q of the RNA hairpin at 
the crossover force / = 17.3pN. 



the value of the attempt frequency fco by fitting (|7I8|) 
to the prediction by Kramers theory Q. We obtain 
fco = 10 5 s _1 which is of the order of magnitude of the 
values reported for other hairpins [lflj . The agreement 
found between the predicted effective barrier 0} and the 
results from the experiments I|7I8|I validates our descrip- 
tion hence providing a way to estimate the attempt fre- 
quency fc 20] . The representation of the effective bar- 
rier as a function of the applied force reflects two dis- 
tinct regimes (Fig|2} characterized by different slopes of 
Bes(f). These correspond to different locations of the 
effective barrier (|5j) and different values of the fragility 
©. We define a crossover force / as the value at which 
the extrapolated straight lines corresponding to regimes 
I and II intersect each other (Fig|2}. 

A kinetic barrier is characterized by its location n* (/) 
and its height G„. (/). As shown in the inset of FigEl 
the free energy landscape at force / shows that there are 
two barriers corresponding to transition states located at 
n* rs 15 - 19 and n\ ps 6 - 7. At low forces, f < f 
(regime I), the highest barrier is located at n* and cor- 
responds to the entropy cost associated to the opening 
of the four bases loop. Whereas for large forces, f > f 
(regime II), the kinetics is governed by the barrier lo- 
cated at n*, at the interface between the GC and AU 
rich regions of the hairpin. In our experiments we ob- 
serve the two different regimes, I and II, because the 
crossover force, / = 17.3pN, is within the experimen- 
tally accessible range of rupture forces, [/ c ± 3pN] (in- 
set of FigflJ, where f c — 17.7pN is the critical force 
verifying Gn(Jc) = Go(/ c ) in Cfjl. In order to inves- 
tigate the unfolding/refolding kinetics over a broader 
range of forces than those accessible in force-ramp experi- 
ments, force-jump experiments |2l| could be very helpful. 
From B^g P and we obtain the fragility fj,(f) by us- 



FIG. 3: Fragility of the RNA hairpin as a function of the ap- 
plied force. The continuous line is the prediction by Kramers 
theory l-SKit whereas symbols correspond to the values ob- 
tained from pulling and relaxing data. We also show the av- 
erage value of the fragility at / c (filled square). Inset: Recon- 
struction of the zero- force free energy landscape G°(n) com- 
pared to the Mfold prediction [Ts|l . Error bars indicate vari- 
ability in the estimates obtained for different pulling speeds. 



ing lloKit . In Fig|3]we show the agreement between the 
fragility obtained from the experimentally measured bar- 
rier H7I8|I and Kramers theory 0}. Finally, our method 
can be used to experimentally reconstruct the free en- 
ergy landscape of the molecule from the sole knowledge 
of the breakage force distribution at different loading 
rates. We first determine the location of the force de- 
pendent transition state n*(f) from the measured value 
of x e g (/) (JSJ). Using the saddle point approximation 
we then identify the effective barrier with the largest 
contribution to the sum appearing in the r.h.s of 0}, 
B°g P (f) pa Max n G„(/) = G„»(/), where we have taken 
h(n*) ps 1 . Finally we apply Q and extrapolate the 
free energy G n *(f) to zero force to obtain G n * (/ = 0). 
In the inset of Fig|3] we show the experimentally recon- 
structed free energy landscape compared with the results 
obtained from Mfold [l5| . 

Changes in the position of the transition state JS} along 
the reaction coordinate axis with the force can corre- 
spond to two different situations: (i) The location of the 
transition state n* does not change with force but the ex- 
tension x n (f) does for all configurations as predicted in 
|23j |: (ii) The free energy landscape of the hairpin shows 
multiple barriers leading to different transition states de- 
pending on the value of the force. This is the case con- 
sidered in the present study. Interestingly the value of 

varies with force in the situation (ii) but not in (i). 
Therefore /i is a good parameter to identify structural 
changes in the transition state. 

A useful analysis of experimental data for molecular 
rupture is the two-states model [2~1 . Ejf where the posi- 
tion of the kinetic barrier along the reaction coordinate 
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is fixed: x F , x VF and the fragility, p, = (x F — x lJF )/(x F + 
x UF ), do not depend on the force. In this approxima- 
tion, ln[r(ln(l/P 11 i y>)) 1 is a straight line as a function of 
the applied force [18|. A plot of the experimental data 
for ln[r(ln(l/P u ( r )))] versus force displays a non-zero cur- 
vature (data not shown) indicating a force dependent 
fragility. Moreover, in two-states systems the dependence 
of the mean value of the unfolding and refolding forces on 
the rate r can be estimated in the experimental regime 

where « 1 M 

</> u oc ^Hr)], </}r oc [ln(r)]. (9) 

By fitting the experimental results to we can estimate 
values for x F ,x VF which give (i « 0. This corresponds 
to a barrier located in the middle between the F and UF 
states, in disagreement with the free energy landscape 
shown in Fig[3 Yet, this value coincides with the average 
fragility measured over the range of forces / G [15 — 
19pN] (FigOJ suggesting that fragility estimates obtained 
by fitting the two-states model to the experimental data 
correspond to averages of force dependent fragilities over 
the range of breakage forces explored in the experiments. 

We have applied Kramers theory to investigate the ki- 
netics of unfolding and refolding of an RNA hairpin under 
mechanical force. The analysis of the experimental data 
for the unfolding and refolding force distributions allows 
us to determine the location of the force-dependent ki- 
netic barrier, the attempt frequency kg of the hairpin and 
the free energy landscape of the molecule. The method 
should be applicable to hairpins with multiple barriers. 
The theory presented here may fail to describe the un- 
folding/refolding of the hairpin at low forces and/or high 
temperatures, where breathing configurations are rele- 
vant [23| and the free-energy landscape becomes mul- 
tidimensional. The presence of force-induced structural 
changes in molecular transition states is a general feature 
of biomoleculcs typically showing a rugged free energy 
landscape. Proper consideration of the force-dependence 
of the fragility is crucial to correctly interpret the results 
from pulling experiments and to relate force unfolding 
measurements with thermal denaturation experiments. 
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